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^_>^ ' Despite the fact that by now one dimensional and three dimensional systems of inter- 

>D , acting particles are reasonably well understood, very little is known on how to go from 

(>_P ■ the one dimensional physics to the three dimensional one. This is in particular true in 

a quasi-one dimensional geometry where the hopping of particles between one dimen- 

C^ ' sional chains or tubes can lead to a dimensional crossover between a Luttinger liquid 

S , and more conventional high dimensional states. Such a situation is relevant to many 

I 1 physical systems. Recently cold atoms in optical traps have provided a unique and con- 

'"^ ' troUable system in which to investigate this physics. We thus analyze a system made of 

r^ ' coupled one dimensional tubes of interacting fermions. We explore the observable con- 

Q , sequences, such as the phase diagram for isolated tubes, and the possibility to realize 

^J ■ unusual superfluid phases in coupled tubes systems. 

Keywords: Cold atoms; Superconductivity; Luttinger liquids. 

X. 

^ . 1. Introduction 

As is well known, interactions have drastic effects on the behavior of bosonic and 
fermionic systems, and change drastically their properties compared to those of 
noninteracting particles. Understanding the properties of such strongly correlated 
systems is a particularly challenging problem. The effects of interaction are also 
considerably enhanced when the dimension of the system is reduced. Very natu- 
rally this quest for strongly correlated systems has thus led to investigations of low 
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dimensional systems, in particular in condensed matter systems. Among these low 
dimensional systems, one dimensional ones play a special role. In such systems, in- 
teraction effects are particularly strong since there is no way for particles to avoid 
each other. For fermions interactions are known to change drastically the properties 
compared to the canonical ones of a Fermi liquid. The one dimensional interacting 
fermionic system is indeed one of the very few solvable case of a non-fcrmi liq- 
uid, known as a Luttinger liquid^. For boson systems, here again the interactions 
lead to properties quite different from the ones of weakly interacting bosons in one 
dimension, and that are much closer to the ones of fermions: the interacting one 
dimensional bosonic systems is again a Luttinger liquid. 

For fermions the progress in material science and nanotechnology have recently 
made possible to realize such one dimensional systems. Another class of materials 
in which such low dimensional systems could be achieved with an unprecedented 
level of control has been recently provided by cold atoms^'"^. Indeed cold atoms not 
only provide the possibility to realize both bosonic and fermionic low dimensional 
systems, but the interactions and kinetic energy can be controlled at will using 
optical lattices and Feshbach resonances^'^'^. 

Such systems allow also to tackle another question that is of importance for a 
large class of materials. Since the effects of interactions vary enormously with the 
dimension, it is crucial to understand how one goes from the one dimensional behav- 
ior to a more conventional three dimensional one. Understanding such dimensional 
crossovers is crucially important for materials such as the organic superconductors^, 
and even more complicated anisotropic materials such as the high Tc superconduc- 
tors. Cold atomic systems thus provide a unique system for which these questions 
could be investigated. This can be done by studying the properties of quasi-one 
dimensional systems made of many one dimensional tubes coupled together. Such 
systems have been both investigated^ and realized'^ for the case of bosons, and it is 
now possible to achieve similar trapping in optical lattices for the fermionic case as 
well^. In these notes we thus focuss on such fermionic systems. We show how the 
peculiarity of cold atoms leads to interesting properties already at the level of an 
isolated ID system. We then examine in details the effect of the coupling between 
the tubes and shows that such systems can be potential realizations of triplet su- 
perconductors. We expand in these notes on Rcf. 10 and refer the reader to this 
paper for technical details and more specialized references. 

2. Fermions, one dimension 

Let us first consider interacting one dimensional fermions. We use the following 
extension of the traditional Hubbard model 

i/ = - ^ tcr (4mC(Tm+l + ll-C.) +U^ n^ranira- (1) 

cr,?n m 

This Hamiltonian describes a ID Fermi gas on a spin dependent periodic potential. 
tc is the hopping from site m to site m + 1 for the spin species a and U the local 
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interaction. In the cold atoms context the spin index refers to two hyperfine states, 
or two different types of atoms {e.g. ®Li and ^"K). Even though there may be no 
true spin symmetry, we use the spin language to describe this binary mixture. It is 
important to point out differences between (1) and the standard Hubbard model. 
First because the system is in an optical lattice it is possible to control separately 
the hopping of the up and down particles. As we discuss below this allows for a 
richer phase diagram than the one of the standard Hubbard model. Second, and 
most importantly, there is no mechanism in the cold atom context that allows to 
relax the "spin" orientation, since this corresponds to two different isospin states. 
As a result the number of spin up and spin down is .separately conserved. This is 
at strong variance with a standard condensed matter situation where one generally 
impose a unique chemical potential for the spin up and spin down particles. The 
important difference, as shown in Fig. 1, is that there are, for the case of equal 
number of spin up and spin down species, only two Fermi points instead of four. 
For such a system, the various operators corresponding to the different types of 




Fig. 1. Relation dispersion E(k) for a system with spin dependent hopping. Left: in the traditional 
case there are spin relaxation phenomena ensuring the equilibrium of the two spin species. There 
is thus a unique chemical potential for both species. This implies the presence of four fermi points. 
Right: for cold atoms, there is no such relaxation and the number of spin up and down are separately 
conserved. For an equal number of spin up and down, there is thus only two Fermi points since 
the two fermi momenta kp-^ = kp ,. 



order in the tubes are of two types^. There are charge density or spin density order 
described by 



Ol(x) ^Y.'^\{x)al^,,^^,(x) 



(2) 
(3) 
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where a^ , a^ , a^ are the PauU matrices The operator Op is simply the total density, 
whereas 0%{x) measures the spin density (along the direction a = x, y, z). We have 
used the notation in the continuum c^ -^ ^{x) as the operator destroying a fermion 
at point X. The density contains a g ^ component and a g ~ 2kF one. The 
operators giving the 2fci? modulation of the charge or spin density are 

OcDw(a;) = i^^R^^L^i^) + Ai^Lii^) = ^-^'''^^" cos(V20„) (4) 

01t,w{^) = A^^Lii^) + Al^L^i^) ^ ^^^^''^^'' cos(V20.) (5) 

OiDw(^) = ~^{A^^Lii^) - V'ljiV'iT(^)) = '^J^'^'^ e'^'^- sin(V20,) (6) 

OiDw(^) = ^I^tV-ltC^) - Ai^Lii^) = ^^^e^^^''2sin(V20^) (7) 

where we have used the notation '0]j (resp. il)]^) to denote the operator creating a 
fermion with a momentum close to +fci? (resp. —kp). These operators have a simple 
representation in terms of continuous field (j)p^^ (resp. Op^a) describing the collective 
excitations of charge and spin densities (resp. current densities). This representation 
is the so-called bosonization representation and is very useful to identify the types 
of order in the system. For more details on this representation see e.g. Ref. 1. Here 
we will only use that (/) and V0 have canonical conjugation commutation relations, 
so order in one of the field implies exponentially decreasing correlations in the other. 
Besides the above mentioned types of order, other instability in the system are the 
superfluid instabilities of the singlet or triplet type 

(7. a' 

where SS denotes singlet pairing whereas TS is triplet pairing. These operators 
describe paring with zero total momentum. Other pairings are of course possible 
but are usually less relevant. These operators become 

Ossix) = i^l^^liix) + ^bl^4'^j^,^{x) ^ -^e-^^^" cos(V20,) (10) 

O^si^) = V'I.tV'ItW +^-[4^1,4(2;) = -le-*^^«p cos(V20,) (11) 

O^s(x) = -ti^l^^pl^ix) ~ i^l^i^^nii^)) - ^e-^^^" sin(V20,) (12) 

Ohi^) = V'ljTV'li(^) - i^h^Rii^) - ^^e-^^^^" sin(V2^,) (13) 

We will not dwell on the solution^" of the model here and will only quote the 
results. Compared to the standard case of the Hubbard model the main effect of 
the spin dependent hopping is to open a spin gap both for repulsive and attractive 
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interactions (the standard Hubbard model being niassless for repulsive interactions) . 
For attractive interactions the spin gap corresponds to the formation of singlets. 
In the bosonization language this corresponds to the order (pu —^ 0. Correlation 
functions containing 6cr or sin(-\/2'/'(T) thus decay exponentially to zero. The leading 
instabilities are thus the charge density wave (CDW) order or a singlet superfluid 
(SS) instability. Physically it means that fermions of opposite spin pair and that 
these pairs behave roughly as bosons and can then condense to give a superfluid 
(SS) or crystallize to give a charge density wave (CDW). 

For repulsive interactions the hopping difference induce an order in 0o- — > vr/vo- 
Correlation functions containing 9a or cos(-\/20cr) decay exponentially to zero while 
sm{^/2(f>cr) tends to a constant. The leading instabilities are thus a spin density 
wave order along the z direction (SDW^) or a triplet superfluid instability along z 
(TSz)- Which one of these instabilities dominates is dictated by the decay of the 
correlations of the charge sector that remains niassless. For repulsive interactions, 
the SDWz one dominates, while the TSz one decays with a larger power law and 
is thus subdominant. A summary of the phase diagram is given in Fig. 2. For 



u 



H 



BKT 



SDW (TS) 



Non-int FG 



SG(?) 



U. 




CDW (SS) 



SS/CDW 



Fig. 2. Schematic phase diagram for the model in Eq. (1) with equal number of spin up and 
down fermions away from half-filling. The interaction strength is U and z = \t^ — t||/{t| + i|). 
All phases (SDW: spin density wave, CDW: charge density wave, SS singlet superfluid, TS triplet 
superfluid) exhibit a spin gap As (however, A^ = for {7 = and z = with (7 > 0. A cartoon 
of the type of order characterizing each phase is also shown. In the area between dashed lines the 
dominant order (either CDW or SS) depends on the lattice filling. In the SG phase, spin up and 
down fermions are segregated (demixed). 
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repulsive interactions the spin sector is thus stiU massive but the spin gap docs not 
correspond to the formation of singlets. On the contrary it is due to a breaking of 
the spin rotation symmetry with formation of an ising-like antiferromagnetic order 
along the z direction. Note that even if the spin sector is massive, the spin density 
wave order is not perfect (except at a commensurate filling of exactly one fermion 
per site) due to the fluctuations in the charge sector. Contrarily to the case of the 
Hubbard model where the SDW order in the three directions decays in a similar 
way (due to the spin rotation symmetry) here the SDWz order decays as a power 
law because of charge fluctuations while SDWa;_y order decays exponentially fast 
because of the presence of the gap in the spin sector. We refer the reader to Ref. 10 
for more details on the various correlation functions and ways to probe the existence 
of such a gap by Raman spectroscopy. 



3. Coupled Fermionic tubes 

Let us now turn to the case of coupled tubes. The coupling between the tubes can 
be described by the single particle hopping term 

H± = -tl_ Y^ Ycl,„i,aC<y,m,f5 (14) 

where a and (3 are the tube indexes and {a, (3) denotes nearest neighbor tubes. 
Because for cold atomic gases the interactions are short range, there is no need to 
take into account interactions between different tubes. If the single particle hopping 
(14) is of the same order of magnitude than the ones in (1) then the system is a 
three dimensional system, while if tj^ = the tubes are uncoupled. The Hamiltonian 
(14) is thus able to describe the dimensional crossover between these two situations. 
For bosons the effects of such a term have been investigated in Ref. 8. For fermions 
treating such a single particle hopping is an extremely challenging problem. Indeed, 
contrarily to the case of boson, the average of a single fermion operator does not 
exist, and thus the term (14) cannot be decoupled simply. Various approximations 
have thus been used to tackle this term and we refer the reader to Ref. 11, 7, 1 for 
details and further references. 

However for the particular case of the fermionic system investigated in these 
notes, one is in a much more favorable situation because of the presence of the 
spin gap. In the case of attractive interactions, the fermions form singlet pairs, that 
essentially behave as bosons. One is thus essentially led back to the bosonic case 
of Ref. 8. We thus focus here on the more interesting repulsive case. Because of 
the presence of the spin gap the situation is now quite different to the one where 
each tube is described by the simple Hubbard model. In the latter case the spin 
excitations are massless and the single particle correlations decay as power laws 

{TrC^,,X^)cl,{0)) ex i^-j (15) 
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where Kp is an interaction dependent parameter (the Luttinger liquid parameter) 
characterizing the isolated tube. In that case a simple scaling analysis to second 
order in perturbation in the single particle hopping term (14) shows that it scales 
as 

i^L4-i[K, + l/K,]-l (^g) 

and is thus in general a relevant perturbation, driving the system away from the 
isolated tube fixed point. However in the presence of a spin gap the single particle 
excitations now decay exponentially 

(T,c„,,, (a;)4,,(0))oce-l-l/« (17) 

where ^ is the correlation length induced by the presence of the spin gap (typically 
of order ^ = vp/^s if vf is the Fermi velocity). The single particle hopping is now 
an irrelevant perturbation to (1). Physically it simply means that because of the 
presence of the spin gap it is now impossible for a single particle to leave a tube since 
it would take away one spin and thus destroy the spin gap. One would therefore 
need a critical value of t± for this to happen. However, although the single particle 
hopping is an irrelevant operator, it can generate relevant perturbations at higher 
order^ . Such relevant perturbations corresponds to hopping between different tubes 
that do not destroy the spin gap. 

Quite generally, as shown in Fig. 3, two types of relevant perturbations are 
possible when looking at the second order term generated by (14). The first one 
corresponds to particle-hole hopping between the tubes, it is thus a density-density 
or spin-spin term, while the second one is the hopping of two particles between 
between the tubes and thus corresponds to a Josephson coupling. The simplest way 
to analyze the terms generated by (14) is to use the bosonization representation. 
The intertube tunnelling term becomes 

o-=±l,r=±l(Q,/3>"' 

Going to second order in (14), using the fact that ip^ is now ordered, and the fact 
that all operators leading to exponential decay must be eliminated to get the leading 
operator, one sees that a first surviving term is 

HsDW = JsDW^ / dx0ljjy^^^^{x)0sDW,j3ix) (19) 

a,f3 •' 

The value of the coupling constant can be obtained by standard second order per- 
turbation theory and is of the order of Jsdw ~ t'j_/^s- This term is thus a standard 
superexchange term coupling the antiferromagnetic spin modulation on two differ- 
ent tubes. Quite naturally this term would tend to stabilize the SDWz instability 
that would develop in an isolated tube and lead to a three dimensional ordered 
antiferromagnetic phase. 
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Fig. 3. Relevant coupling generated to second order by the single particle intertube hopping. In 
order to avoid breaking the spin gap, a correlated hop of two objects must be done simultaneously. 
Left: particle-hole hopping. Since this term exchanges two particles on two different tubes, it 
corresponds to a standard density-density or spin-spin coupling term. This is the superexchange 
term. Right: particle-particle hopping. This term allows a pair to jump from one tube to the next 
and is thus analogous to the Josephson coupling between superconductors. 



There is however another term that survives. This is a term where two particles 
can hop from one tube to the other. It is of the form 



Hts 






dxO\,g^ {x)Ots,A^) 



(20) 



This term corresponds to the hopping of a pair of fermions, that are in a triplet 
paring state from one tube to the other. Singlet hopping is here cancelled because 
of the presence of the spin gap. Because a superfluid pair has a global momentum of 
zero, the pair hopping does not contain the oscillating 2kp factor that is present in 
the superexchange term (see (4)). Here again the coupling constant is of the order 
oiJTs-tl/As. 

Both the superexchange and the Josephson term tend to stabilize their corre- 
sponding type of order. What phase is realized normally depends crucially on what 
is the leading one dimensional instability. Usually the Fermi momenta of all tubes 
are the same and there is no oscillatory factors in (19). Since the SDW order is 
dominant for the isolated tube, and both Jsdw a-nd Jxs are of the same order of 
magnitude a simple RPA treatment of the coupling term shows that the SDW phase 
is stabilized. This is the situation depicted in Fig. 4(a). This situation is normally 
the one realized in condensed matter systems. It corresponds to the stabilization 
of a three dimensional antiferromagnetic order for repulsive interactions. However 
the situation can be much richer if the Fermi momenta of the different tubes are 
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different. Note that for cold atoms this is the rule rather than the exception be- 
cause of the presence of the parabolic confining potential that makes each tube 
different depending from its distance from the center. This situation can also be 
reinforced artificially by adding an additional modulation of the optical lattice. In 



kp^O which means that the oscillatory factors remain 



that case one has kp 
in (19). Such factors considerably weaken the intertube superexchange. Physically 
this means that the antiferromagnetic fluctuations on the neighboring tubes are now 
incommensurate with each other and thus cannot couple very well as depicted in 
Fig. 4(b). On the other hand the Josephson term which is a q ~^ transfer between 
the tube is not affected by such a difference of Fermi momentum and remains un- 
changed, as schematically shown in Fig. 4(c). One is now in a situation where this 
term can become dominant^" even if the triplet superfluid instability in the isolated 
tube is a subdominant instability. One would thus be in a situation to stabilize a 
three dimensional triplet superfluid phase. 






(a) 



' Q d Q ^ 



(b) 



CM300" 




(c) 



Fig. 4. Schematic diagrams for the true long range order induced by inter-tube particle-pair or 
particle-hole hopping, a) Antiferromagnetic order for f7 > and different velocities, with equal 
Fermi momentum (i.e. equal number of fermions in each tube), b) Frustration (incommensurabil- 
ity) to the antiferromagnetic order between tubes when the Fermi momentum is different between 
the tubes, c) No such frustration effects when the order parameter carries zero total momentum, 
namely for the superfluidity. 
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This is a rather unique situation. From the theoretical point of view such a 
stabihzation of a subdominant instabiUty by intertube coupUng is quite interesting 
and potentiaUy relevant to other situations as well. In particular, this mechanism is 
similar, albeit for a triplet superconductor, to the one advocated to stabilize singlet 
superconductivity for coupled stripes^^. Such cold atomics system could thus be 
controlled systems in which to check for the feasibility of such a mechanism. More 
directly it would be very interesting to have a realization for a triplet supercon- 
ductor. Triplet superconductivity is indeed quite rare. Besides Helium 3, strontium 
ruthenate is the only one candidate well identified in condensed matter^^. Organic 
superconductors, another system made of coupled chains, is also a potential candi- 
date for such unusual superconductivity^'*. Cold atomic gases of fermions could thus 
help shed a light on the mechanisms and properties of unusual superconductivity 
in these very anisotropic systems. 
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